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ABSTRACT 

The recent discovery that the close-in extrasolar giant planet, IID209458b, transits its star has provided 
a first-of-its-kind measurement of the planet's radius and mass. In addition, there is a provocative 
detection of the light reflected off of the giant planet, r Boo b. Including the effects of stellar irradiation, 
we estimate the general behavior of radius/age trajectories for such planets and interpret the large 
measured radii of HD209458b and t Boo b in that context. We find that HD209458b must be a hydrogen- 
rich gas giant. Furthermore, the large radius of a close-in gas giant is not due to the thermal expansion 
of its atmosphere, but to the high residual entropy that remains throughout its bulk by dint of its early 
proximity to a luminous primary. The large stellar flux does not inflate the planet, but retards its 
otherwise inexorable contraction from a more extended configuration at birth. This implies either that 
such a planet was formed near its current orbital distance or that it migrated in from larger distances 
(>0.5 A.U.), no later than a few times 10'' years of birth. 

Subject headings: stars: individual (IID209458, t Bootis) — (stars:) planetary systems — planets and 
satellites: general 



1. INTRODUCTION 

The recent indirect detections of extrasolar giant plan- 
ets (EGPs) by Doppler spectroscopy have taught us that 
such planetary systems can be very much unlike our own. 
To date, nearly '- "30 EGPo have been diocovcrc d around ■ 
Jtaro with pp ootml typop from H i ll to F7 (e.g., [P i layor 4 " 



Queloz 1995 



Marcy et al 19994 pscher et al. 1999| ). The 



than that of Jupiter itself (Guillot et al. 1996). This pre- 
diction has recently been vcrihed by the observation of 
the transit by HD209458b of its primary. The depth of 
the transit yields values for its rad ius that rang e from 
^1.27 Rj (iCharbonncau et al 1999b| ) to -1.7 Rj ([Henry 
et al. 1999| ), with a best value near —1.4 Rj ( Mazel 
et al. 20001). Importantly, since HD209458b transits its 



be- 



planets themselves have minimum masses (MpSin(i ) ) 
tween —0.42 Mj and —10 Mj , orbital semi-major axes 
from —0.042 A.U. to —3.8 A.U., and eccentricities as high 
as —0.71. Such variety vastly expands the parameter space 
within which both theorists and observers must operate in 
defining the physical character of extrasolar planetary sys- 
tems. 

The most interesting, unexpected, and problematic sub- 
class of EGPs are those found within —0.1 A.U. of their 
primaries, 50-100 times closer than Jupiter is to our Sun. 
At such orbital distances, due to prodigious stellar irra- 
diation alone, an EGP can have an effective temperature 
(Teff ) greater than 600 K. Indeed, the EGPs HD187123b, 
HD209458b, r Boo b, HD75289b, 51 Peg b, v And b, and 
HD217107b hkely aU have Tcff s above 1000 K. This is to be 
compared with ToftS for Jupiter and Saturn of 125 K and 
95 K, respectively. The minimum masses of these close-in 
EGPs imply that they occupy the high-Tcff /low-gravity 
corner of parameter space for which t he compositions and 
composition profiles may be un ique (Seager and Sasselov 
19981 ; iGoukcnleuque al 1999| ). 

However, high stellar fluxes on a close-in EGP can have 
profound structural consequences for the planet. In par- 
ticular, stellar insolation can be responsible for maintain- 
ing the planet's radius at a value 20% to 80% larger 



primary, astronomers can derive sin(z), from which the 
planet's mass can be directly determined. This is a major 
advance in the emerging study of extrasolar planets. In 
this paper, we focus on the HD209458 system and what 
broadly can be concluded theoretically from these new 
transit data. We do not provide detailed models from 
which one can extract bulk or atmospheric composition. 
Rather, we show that the large radius of HD209458b is 
a consequence of the retardation by stellar irradiation of 
the otherwise natural cooling of the convective core of the 
planet and that such a radius requires that the planet did 
not dwell for long, if at all, at large orbital distances after 
its formation. Hence, the large radius of a close-in EGP is 
not due to the thermal expansion of its atmosphere, but 
to the high residual entropy that remains throughout its 
bulk as a consequence of its early proximity to a luminous 
primary. 

Recently, Cameron et al. (1999), using spectral decon- 
volution, claim to have seen r Boo b in reflection. To 
investigate this exciting possibility, we calculate a range 
of theoretical radii for the r Bootis planet. As we show, 
the Cameron et al. value of 1.6-1.8 Rj for r Boo b's radius, 
if verified, is a challenge to the still embryonic theory of 
massive, close-in EGPs. Nevertheless, with the discovery 
of a transiting planet, HD209458b, with the maturation 
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of the technique of spectral deconvolution, and with the 
anticipated development of adaptive optics and in tcrfer- 
ometry for the direct study of extrasolar planets (Angej 
1994|] we are clearly entering a new phase in extrasolar 
planetary research. 

2. A SUMMARY OF THE DATA 

The transit of the F8V/G0V star HD209458 (at a dis- 
tance of 47 parsecs) by HD209458b lasts ~3 hours (out of 
a total period of 3.524 days) and has a depth of ^1.5 -2.0% . 
Its ingress and egress phases each last ^2 5 minutes (Char 



bonncau et al. 1999b; Henry et al. 1999). The properties 



of the planet, in particular its orbital distance and its ra- 
dius, scale with the properties of the star and the most 
recent study of HD209458 was conducted in the context of 
these transits by Mazeh et al. (2000). They conclude from 
/ogg/Tcff spectral-line fits and Mv/{B — V) photometric 
fits that l.liO.l Mq, = 1.2±0.1 Rq, Toh ~6000 
K, [Fc/H]-0.0, t = 5.5±1.5 Gyr, and L, ^ 2.0 L© . Mazeh 
et al. then derive for the planet: Rp — 1.40±0.17 Rj , Mp= 

0.69 ± 0.05 Mj (cxM^/^), i = 86°. 1 ± 1°.6, and a = 0.047 

1 /3 

A.U. (ocM* ). (All symbols have their standard mean- 
ings.) 

Fuhrmann et al. (1998) provide parameters for the F7V 
star T Bootis: M*= 1.42 ± 0.05 Mq, R* = 1.48 ± 0.05 
Rq, Teff -6360 K, [Fe/H]~-h0.27±0.08, i = 1.0 ± 0.6 Gyr, 
and L» 3.2 L© . Its Hipparcos distance is ^^15. 6 par- 



secs. iWlth Llic Fuhi'maim tt a/, mays foi' t Bootis, But- 
ler et ' al. (1997) would have obtained foi its close-in EGF. 
MpSin(i) = 4.33 Mj , a - 0.049 A.U., and P ^ 3.313 
days. Charbonneau et al. (1999a) quote an upper limit 
at '^4900A of 5 X 10~^ to the fraction of the star's light 
reflected off the planet. Cameron et al. (1999) claim to 
have detected in the blue-green region of the spectrum 
a reflected fraction of 1.9 ± 0.4 x 10~^. From the semi- 
amplitude ('^74 km s^^) of the Dopplcr shift of this re- 
flected fraction, they derive an orbital inclination (?) of 
29°, which would yield a mass for r Boo b of ~9 Mj . 
From their reflected fraction, an orbital planetary phase 
function, and a geometric albedo {Ag) of 0.55 (similar to 
that of Jupiter in the visible), Cameron et al. obtain a ra- 
dius for r Boo b of 1.6-1.8 Rj . If Ag were smaller, the 

inferred radius of the planet would be larger (oc Ag ). 

3. THE RADH OF CLOSE-IN EGPS 

By whatever processes giant planets are initially assem- 
bled, they must start out significantly larger than they end 
up. In isolation, they would cool inexorably due to radi- 
ation from their surfaces and shrink accordingly, just as 
does a protostar or a brown dwarf. Early on, due to the 
negative effective specific heat of an object in hydrostatic 
equilibrium supported by ideal gas pressure, energy loss 
results in an increase in its central temperature. However, 
the density increases more quickly and, as a consequence, 
the specific entropy [S) monotonically decreases. Since 
EGPs are almost f ully convective (ev en if under significant 
stellar irradiation; puillot et al. 1996 ), they are isentropes. 
This is an essential point. Given an EOS and a planetary 
mass, an EGP's core entropy determines its radius (and 
its surface gravity). A large radius is a consequence of a 
large entropy. 

As Zapolsky and Salpeter (1969) demonstrated for plan- 



ets made of high-Z material, any planet with a cold ra- 
dius larger than —0.5 Rj , must be made predominantly 
of hydrogen. Usi ng the ANEOS equation of state tables 
(Thompson 1990), we derive that an "olivine" (rock) or 
H2O (ice) planet with a mass of 0.69 Mj has a radius of 
0.31 Rj or 0.45 Rj , respectively. Importantly, these radii 
are 3-4 times smaller than observed for HD209458b and 
prove that HD209458b must be a hydrogen-rich gas giant; 
it cannot be a giant terrestrial planet or an ice giant such 
as Neptune or Uranus. 

The rate with which an EGP shrinks is determined by 
the opacities in its outer radiative zone and the degree of 
stellar irradiation. If there were no radiative losses, the 
EGP would not shrink. In isolation, the energy loss rate 
and Teff are determined in the context of a self-consistent 
radiative/convective model and Tes is a function of only 
S and Mp. (Note that, for a given metallicity, the sur- 
face gravity of the EGP is a function of S and Mp alone.) 
The atmospheric flux in the outer skin determines the 
temperature profile at the boundary of the convective 
core and, hence, the rate of c ore entropy and radi us de- 
crease. According to theory ( |Saumon et al. 1996| ; Bur- 
rows et 1995,1997), isolated EGPs shrink rapidly. A 
1-Mj EGP in isolation contracts below 2.0 Rj in less than 
10^ years. The theory depends upon the atmospheric 
opacities and metallicity an d can reproduce the current 
Jupiter (Hubbard et al. 1999), but there still remain m ajor 



uncertainties co ncerning grain and cloud formation (Lu 
nine et al. 198S| ) , rainout ( Burrows and Sharp 1999 ) , gas 
phase opacities ( Burrows et al. 19971), and the depth of 



penetration of the stellar flux ( jCuillot et al. 2000 ) . Hence, 
while the basic theory of EGP evolution is firm, the de- 
tails are not and significant ambiguities in the variation of 
Rp and Tes with age persist. 

It is by altering the temperature/pressure profile of the 
atmosphere of an EGP that stellar irradiation can retard 
the evolution of the EGP's core entropy and, hence, Rp . 
Essentially, irradiation flattens the temperature profile at 
the top of the convective zone, while at the same time 
moving the radiative/convective boundary inward. The 
consequent growth of the radiative z one is a central fea - 
ture of the large-radius phenomenon (Guillot et al. 1996). 
Since radiative fluxes are governed by the product of ther- 
mal diffusivities and temperature gradients and since the 
thermal diffusivity decreases with increasing pressure, the 
flux of energy out of the convective core and the rate 
of core entropy change are reduced. For close-in EGPs, 
Toft stabilizes early at large values, while the "effective" 
Toff of the core, where most of the heat and mass resides, 
is drastically lowered. The upshot is a retardation of the 
contraction of the planet. Large EGP radii are a conse- 
quence of such retardation, and not of the expansion of the 
outer envelope by stellar heating. This is easy to demon- 
strate by noting that the scale height of HD209458b's at- 
mosphere, under the assumption that Teff is between 1200 
K and 1700 K, is only —1% of Rp . Though the mapping 
between S and Rp is unaltered by irradiation, the map- 
ping between S and age can change significantly. Guillot 
et al. (1996) predicted this behavior for 51 Peg b, using 
a Bond albedo {Ab) of 0.35, obtaining Rp 's of 1.1-1.3 
Rj after 8 Gyr. Scahng the results of that paper using 
the stellar flux on HD209458b and a mass of 0.69 Mj , we 
obtain Rp 's between 1.4 Rj and 1.6 Rj for ages between 
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10 and 3 Gyr, perfectly in line with the transit observa- 
tions (§2). Figure 1 depicts two possible evolutionary tra- 
jectories for HD209458b, if born and fixed at 0.049 A.U. 
Also included on Figure 1 is a theoretical Rp -t trajectory 
for a 0.69-Mj EGP in isolation. We have used for these 



calcul itions the formalisms of Guillot et al. (1995, 1996), 
Guillot and Morel (1995), and Burrows et al. (1997). Tidal 
and Roche effects have been ignored and are relevant only 
at very early ages (^10^ yrs). Note, however, that despite 
its small orbital distance, HD209458b is beyond the Roche 
limit by a factor of ~2.5 and is stable against loss both by 
thermal (Jeans) escape and non-t hermal proces ses involv- 
ing absorption of stellar UV flux ( Trilling 1998| ) . A box of 
empirical ages and Rp 's (Mazeh et al. 2000) is superposed. 

As can be seen in Figure 1, there is a great differ- 
ence between the theoretical radius of an isolated and 
an irradiated EGP. Importantly, since the scale-height ef- 
fect is miniscule, the Rp -t trajectory of the isolated EGP 
(model I) immediately suggests that if HD209458b were 
allowed to dwell at large orbital distances (>0.5 A.U.) for 
more than a few x 10^ years, its observed radius could 
not be reproduced. As Figure 1 demonstrates, it is at 
such ages that the radius of an isolated 0.69-Mj EGP falls 
below HD209458b's observed radius. Note that for >0.5 
A.U. the intrinsic luminosity of our isolated HD209458b 
model at such an age exceeds the amount of stellar light 
that would have been intercepted. Hence, stellar irradi- 
ation would have had little effect on this model. How- 
ever, many orbital-distance/age histories can be co ntem- 



plate d and the se will be the subject of a future work (Guil 
lot et\al. 2000D. To reiterate, HD209458b could not have 



cooled off and achieved a radius or an entropy like that of 
Jupiter and then moved in. If it migrated, it had to have 
moved in early in its life. The radius of a close-in EGP de- 
pends upon its history (and the history of the star); large 
radii require early proximity to a central star. This fact 
provides an upper limit to the timesc ale of planetary m i- 
gration, if migration did indeed occur ( Trilling et al. 1998| ): 
conservatively, HD209458b dwelled less than a fewxlO^ 
years at more than 0.5 A.U. 

The actual evolution of a close-in EGP's radius depends 
upon its B ond albedo (Mar ley et al. 1999 ; Seager and Sas-| 
selov L998; Budarsky, Burrows, and Pinto 1999), the level 
of any clouds formed and their optical depths, the gas- 
phase abundances and opaciti es, the helium and metallic- 
ity fr actions, the H-He EOS (Saumon, Chabrier, fc Van| 
Horn 199E ; ^aumon et a/.199S ), the depth of stella r flux 



penet ration, a nd the primary star. In a later paper (Guil 
lot et al. 2000| ) , we explore these effects and conduct a de- 



tailed parameter study. A range of trajectories similar to 
those depicted in Figure 1 are still possible and when data 
on the HD209458b transit and primary star are further re- 
fined, theorists may well be able to sharply constrain the 
character of HD209458b's atmosphere and composition. 

3.1. Theoretical Radius for t Boo b 

We provide in Figure 1 two representative theoretical 
Rp -t trajectories for r Boo b, along with an error box 
constructed using Cameron et al. (1999) and Fuhrmann 
et aZ. (1998). Depending upon L* and its Bond albedo, 
which might vary from ~0.0 to perhaps 0.6, r Boo b's 
Toff is between 1350 K and 1750 K. These Toffsare higher 
than those corresponding to HD209458b and reflect t 



Bootis' higher L*. Due to uncertainties in the relative 
position of the silicate and iron cloud decks and the re- 
gion of neutral alkali metals, there are still uncertainties 
in the atmospheric composition and albedos of such an 
EGP dMarlcy et al. 1999|; [Sudarsky, Burrows, and Pintc 



1999). A high Bond albedo might be a consequence of the 



presence at altitude of reflecting clouds. Without these 
clouds, both the Bond albedo and the geometric albedo in 
the blue-green region of the planetary spectrum would be 
low, perhaps below 0.1. A low geometric albedo would put 
the measured radius of r Boo b even higher than the range 
quoted by Cameron et al. . As the evolutionary models on 
Figure 1 suggest, though a lower Bond albedo results in 
a larger Rp at a given age, the Rp -t error box for t Boo 
b given in Figure 1 still seems out of reach. Note that a 
9-Mj T Boo b model in isolation achieves a radius of 1.6 
Rj within a scant 5 Myr and that at 1 Gyr such a model 
has a radius of 1.1 Rj . 

The early radii of r Boo b shown on Figure 1 are so much 
smaller than those of HD209458b because of the larger in- 
ferred mass of r Boo b and the particular choice of opac- 
ity data used for these exploratory models. Other opacity 
databases would give quantitatively different radius-age 
trajectories early on, but would not alter the conclusion 
that the theoretical radius of a massive planet at late times 
(> fewxlO'' years) is significantly below the r Boo b error 
box shown in Figure 1. Furthermore, under the assump- 
tions that the Bond albedo is zero, that there is no outer 
radiative zone, and that the atmosphere and core of the 
planet have the same entropy, we can obtain a strong up- 
per bound to the radius of a planet of a given age and mass 
(Guillot et al. 1996). For r Boo b, this radius is 1.58 Rj for 
7 Mj and 1.48 Rj for 10 Mj , both below the quoted ra- 
dius range. Hence, we have difficulty fitting the Cameron 
et al. (1999) refiection data. In particular, despite signif- 
icant irradiation, the large planetary mass measured by 
Cameron et al. (7-10 Mj ) results in rapid early contrac- 
tion. If the Bond albedo is lower, the theoretical r Boo 
b radius for an age near 1 Gyr would be larger, but the 
measured radius would also be larger. 

4. THE EFFECTS OF THE EQUATION OF STATE 

Our cooling calculations use the hydrogen/helium EOS 
of Saumon, Chabrier, and Van Horn (1995, SCVH) . Recent 



shock-compression expe ri ments on deuterium ( Holmes 
Ross, and Nellis 1995|; pa Silva et a ?. 1997] ; |Coll 



et al. 1998 ) show that this EOS underestimates the de- 
gree of molecular dissociation for pressures in the range 
0.1 ^ P ^ 2Mbar. However, by softening the repul- 
sive part of the H2-H and H-H potentials, the EOS of 
SCVH can be made to reproduce all experimental results 
( Saumon et al. 1999| ). Since adiabats of hot Jupiters com- 
puted with the modified EOS are systematically cooler and 
denser than those of SCVH by up to 6% in T and 10% in p, 
such changes in the EOS will reduce the theoretical radius 
of HD209458b by only a few percent, while modifying the 
corresponding quantity for r Boo b by a negligible amount. 
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5. ATMOSPHERIC TRANSMISSION AND REFRACTION 



6. CONCLUSIONS 



Details of the transit lightcurve depend upon the distri- 
bution of slant optical depth, as determined by molecular 
opacity and cloud layers, and, to a lesser extent, on refrac- 
tive redistribution of the stellar surface brightness by the 
planetary atmosphere. 

We use a theory to compute refractive effects that is 
essentially identical to that of the standard the ory for oc- 
cultations of stars by p lanetary atmospheres ( Hubbard] 



Yelle, and Luninc 1990) and have calculated the stellar 



brightness distribution for the transit of a planet with a 
radius equal to 1.4 Rj and a hot solar-composition atmo- 
sphere. The slant optical depth r is computed for molec- 
ular opacity sources (with and without clouds) and at a 
variety of wavelengths. The effective radius of the planet, 
as determined by fitting the transit lightcurve, will de- 
pend upon the gas-phase opacity, the molecular composi- 
tion, and the location of t he optically-thick cloud layers 
( Seager and Sassclov 199"9| ). It will also be a diagnostic 
function of wavelength. Our preliminary calculations for 
cloud- free models indicate that |Ai?p/i?p| between 4500 A 
and 6500 A might be as much as 3%, smaller (~G.5%) if the 
alkali metals are i mportant. Using a s i mple model for sil i- 



cate cloud growth (Lunine et al. 1989 
we find a cloud base for HD209458 



Marley et al. 1999] ), 
between 10^^ and 



0.4 bars and grain particle sizes between 1 and 100 mi- 
crons. While for cloud models the wavelength dependence 
is muted, since the opacity varies strongly with particle 
size, there exists the possibility of remote sensing of cloud 
properti es with high-precision measurements of upcoming 
transits ( Hubbard et al. 2000| ). 



Our theoretical calculations allow us to draw several spe- 
cific conclusions: 

1. HD209458b is a real object, made predominantly of 
hydrogen. 

2. HD209458b's radius is a consequence of the retarda- 
tion of contraction by stellar irradiation and is not 
due to atmospheric expansion by stellar heating. 

3. A large radius such as that of HD209458b requires 
early proximity to its central star. 

Curiously, given L,, a, and Rp , HD209458b's total lu- 
minosity is ^ 1.5 X 10"'* Ijq , about twice that of a star 
with '^lOO times the mass at the very bottom of the stel- 
lar main sequence. Given the large inferred mass of r 
Boo b, its large radius is less easy to explain theoretically. 
However, the inherent difficulties of close-in EGP model- 
ing may yet be responsible for theoretical surprises of a 
qualitative nature. 
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4214, NAG5-7211, NAG2-6007, NAG5-4988, and NAG5- 
4987, as well as by an NSF CAREER grant (AST-9624878) 
to M.S.M. 
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Fig. 1. — Theoretical evolution of the radii of HD209458b and r Boo b (in Rj ; 7 X lO* km) with age (in Gyrs). Model I (solid) is for a 
0.69 Mj object in isolation (Burrows et al. 1997). Models A {Ab = 0.0; Tcff --leOO K) and B {As = 0.5; T<,ff ~1200 K)) are for a close-in, 
irradiated HD209458b at its current orbital distance from birth, using the opacities of Alexander and Ferguson (1994). Models C (Ag = 0.0; 
Teg ~ 1750 K; Mp= 7 Mj ) and D {Ag = 0.5; T(;fr'~1350 K; Mp= 10 Mj ) are for a close-in, irradiated r Boo b, using a similar opacity 
set. The formalism of Guillot et al. (1996) was employed for models A-D. The ranges spanned by models A and B and by models C and D 
for HD209458b and t Boo b, respectively, reflect the current ambiguities in the observations and in the theoretical predictions due to cloud, 
opacity, composition, and depth of flux penetration uncertainties. Superposed are an error box for HD209458b (far right) using the data of 
Mazeh et al. (2000) and one for t Boo b using the data of Cameron et al. (1999) and Fuhrmann et al. (1998). 



